Fish biodiversity declines with dam development in the Lower Mekong Basin

Hydropower dams are a source of renewable energy, but dam development and hydropower generation negatively affect freshwater ecosystems, biodiversity, and food security. We assess the effects of hydropower dam development on spatial–temporal changes in fish biodiversity from 2007 to 2014 in the Sekong, Sesan, and Srepok Basins—major tributaries to the Mekong River. By analyzing a 7-year fish monitoring dataset, and regressing fish abundance and biodiversity trends against cumulative number of upstream dams, we found that hydropower dams reduced fish biodiversity, including migratory, IUCN threatened and indicator species in the Sesan and Srepok Basins where most dams have been constructed. Meanwhile, fish biodiversity increased in the Sekong, the basin with the fewest dams. Fish fauna in the Sesan and Srepok Basins decreased from 60 and 29 species in 2007 to 42 and 25 species in 2014, respectively; while they increased from 33 in 2007 to 56 species in 2014 in the Sekong Basin. This is one of the first empirical studies to show reduced diversity following dam construction and fragmentation, and increased diversity in less regulated rivers in the Mekong River. Our results underscore the importance of the Sekong Basin to fish biodiversity and highlight the likely significance of all remaining free-flowing sections of the Lower Mekong Basin, including the Sekong, Cambodian Mekong, and Tonle Sap Rivers to migratory and threatened fish species. To preserve biodiversity, developing alternative renewable sources of energy or re-operating existing dams to increase power generation are recommended over constructing new hydropower dams.

additional hydropower dams are under construction, including three on the mainstem Mekong River in Laos. Another 19 large dams (> 200 MW) have been proposed across the LMB. By 2040, estimated economic gain from all projects in the LMB will exceed $160 billion USD, while their potential cost is about $170 billion USD due to declining fisheries, forests, wetlands and mangroves 14 . Overall, a net economic loss is anticipated from hydropower projects, with the vast majority of benefits accruing to foreign companies and the majority of losses borne by residents of Cambodia, Thailand, and Vietnam, especially rural subsistence fishers who depend on fisheries for food, protein, and livelihoods 10 .
The impact of dams on the Mekong's fish and fisheries have been investigated and summarized in several recent studies. Ziv et al. 7 estimated fish biomass and biodiversity losses based on future dam alternatives and estimated trade-offs between power production and impacts to fish biodiversity. Kano et al. 15 modeled the future distribution of fish with intense hydropower development and found reduced species richness and available habitat. Ngor et al. 9 assessed the impacts of flow alteration on fish diversity, and reported reductions in local fish richness and abundance in the Sekong, Sesan and Srepok Basins (3S Basin) and altered flow seasonality in the 3S-Mekong-Tonle Sap system, with distinct variations in fish assemblages 16 . Yoshida et al. 12 reviewed the impacts of mainstem hydropower dams on agriculture and fisheries resources, and found that annual fish catch has declined by 276,847 and 178,169 tonnes in Cambodia and Vietnam, respectively. Baird and Green 17 and Green and Baird 18 reviewed and discussed the impacts of dams on fish and ecosystems and provided options for future clean development. Soukhaphon et al. 19 studied several Mekong hydropower dams and found that social and environmental impacts are often cumulative and occur over large scales. All of these studies highlight the significant and varied costs from hydropower development in the LMB and together suggest that further dam construction is not a sustainable development trajectory.
One of the hotspots for hydropower dam development is the Sekong, Sesan and Srepok Basins (3S Basin), which together form a major tributary to the Mekong. In 2007, the Sekong Basin had two dams, the Srepok Basin had two dams, and the Sesan Basin had six dams. In 2014, the number increased to six dams in the Sekong Basin, 13 dams in the Srepok Basin, and 19 dams in the Sesan Basin. Of those dams, six were major dams located on the mainstem river of the Sesan Basin, five were located on the mainstem of Srepok Basin, and one dam was located on the mainstem river of the Sekong Basin-although it was situated far upstream (Fig. 1a). By 2021, 51 dams were operational in the 3S Basin (Fig. 1b), with total generating capacity of 4684 MW. Five more dams are under construction and 79 are planned 13 . Despite the widespread impacts of hydropower dam development, limited research has quantified spatial and temporal changes in fish biodiversity with dam development.
The 3S Basin has vast ecological importance, with a seasonal flood-pulse that fuels exceptional aquatic biodiversity. Moreover, the 3S Basin is one of the main fish migration routes in the LMB, providing spawning and rearing habitats for some species 9,20 . Among the three Basins, the Srepok was historically utilized by 81 species of migratory fish, the Sekong by 64 species, and 54 species of migratory fish used the Sesan Basin 21,22 . The 3S-Mekong-Tonle Sap is a migratory corridor linking fish spawning and rearing habitats, including endangered species like the Critically Endangered Mekong Giant Catfish (Pangasianodon gigas) and Endangered Mekong River Catfish (Pangasianodon hypophthalmus). While the Mekong is home to at least 150 endemic aquatic fauna 23 , 103 of which are endemic fish species living in the LMB 24 , it also has approximately 20 non-native species 25 , which compete with native species for resources 3,24 .
Here, we assess spatial-temporal changes in fish biodiversity from hydropower dam development in the 3S Basin from 2007 to 2014. By using long-term monitoring data, we (1) quantify changes in fish biodiversity composition and indicator species in the 3S Basin, (2) investigate the temporal trends in fish biodiversity (richness, abundance and diversity indices of total, migratory and threatened fish species) and (3) quantify the impacts of dam development on fish biodiversity. We conclude with discussion about sustainable development futures for the 3S Basin.

Results
Overall fish biodiversity and spatial-seasonal changes. Our study area had 247 fish species belonging to 43 families and 14 orders. The most common order in the 3S Basin was Cypriniformes (52%, 129 species), followed by Siluriformes (24%, 60 species), and Perciformes (11%, 28 species). The remaining 11 orders made up less than 3% of the total species counts. PCA analysis demonstrated distinct fish communities among the 3S rivers (Fig. 2a) and changes in those fish communities between the wet and dry seasons in the Srepok and Sesan Basins (Fig. 2b).
We also characterized spatial-seasonal communities in each river by indicator species (Table 1). As whole, Sekong Basin had 19 indicator species, five of which were migratory and two were threatened species (Mystus bocourti-Vulnerable and Pangasianodon hypophthalmus-Endangered), whereas Srepok Basin had 10 indicator species, including 6 migratory and one Critically Endangered species (Probarbus jullieni). Sesan Basin had only one indicator species (Acantopsis sp.). Fish communities in Sekong Basin were relatively similar between seasons, as indicated by the PCA ordination map (Fig. 2b). In the Srepok Basin, there was one indicator species for the dry season and 6 indicator species in wet season. The Sesan Basin had 5 indicator species that characterized fish communities in the dry season and none in the wet season ( Table 1).
The monthly richness of the three categories was not significantly different between the Sekong (42.2) and Srepok (45.2) Basins, but monthly richness was lower in the Sesan Basin (24.1) (Fig. 3b, Table 2). Monthly abundance was significantly higher in the Srepok Basin for migratory fish and in the Sekong Basin for threatened fish (Fig. 3c, Table 2). For the three fish categories, diversity were different among the basins, and diversity was highest in the Sekong Basin, followed by the Srepok Basin, and was lowest in the Sesan Basin (Fig. 3d).
In the 3S Basin as a whole, 8 LMB-endemic fish species were observed, with seven, eight, and four species observed in the Sekong, Srepok and Sesan Basins, respectively. Four non-native species were detected throughout the 3S Basin. The Sekong and Srepok Basins each had three non-native species, while there were none in the Sesan Basin over the study period. Refer to Table S1 for detailed taxonomy.
Temporal fish biodiversity trends and the impact of dam development. Overall, fish biodiversity metrics are increasing in the Sekong Basin and declining in the Srepok and Sesan Basins. Linear regression models indicated that fish biodiversity metrics (species richness, abundance, and diversity) in the Sekong Basin significantly increased from 2007 to 2014 for total fish and migratory fish. Threatened fish biodiversity metrics showed an upward trend for the Sekong Basin, but only diversity was significant (Fig. 4, Supplemental Information Fig. S1). In the Srepok Basin, fish richness, abundance, and diversity significantly declined from 2007 to 2014 for threatened species, while richness also significantly declined for total and migratory species. www.nature.com/scientificreports/ In the Sesan Basin, richness and diversity significantly declined for total species and only diversity significantly declined for migratory species (Fig. 4, Supplemental Information Fig. S1). When the change in biodiversity metrics were significant, model performance (R 2 ) ranged from 0.06 to 0.42 ( Fig. 4, Supplemental Information Fig. S1). These significant but low R 2 values mean that there was considerable noise in the data, or in other words, fish catch varied considerably among months and years and the trend for fish diversity was statistically significant. The highest R 2 values indicated increasing total and migratory fish biodiversity in the Sekong Basin-where few dams were developed, and declining biodiversity for threatened fish in the Srepok Basin.
Based on linear regression at an annual timestep, fish biodiversity metrics in the Sekong Basin increased significantly through time as dams were built in other tributaries of the 3S Basin, primarily the Sesan and Srepok Basins (Figs. 1, 5). A similar pattern was found when the data was partitioned into seasons, particularly in the dry season (Figs. S2, S3). Total fish richness in the Sekong Basin significantly increased from 33 species in 2007 to 56 species in 2014 with a model performance of R 2 = 0.75 (P < 0.01). Abundance, richness, and diversity always significantly increased in the Sekong Basin for total species and migratory species, with R 2 values ranging from 0.47 to 0.79. In contrast, biodiversity trends declined in the Srepok and Sesan Basins as more dams were built (Fig. 5). Monthly total fish richness in the Srepok and Sesan Basins decreased from 60 to 42 species and 29 to 25 species, respectively, from 2007 to 2014. In the Srepok Basin, richness and diversity of threatened fish species significantly decreased (R 2 = 0.56 and R 2 = 0.82, respectively); while in Sesan Basin, total fish abundance and diversity, and the migratory fish diversity significantly decreased (R 2 = 0.51, R 2 = 0.69 and R 2 = 0.63, respectively) (Fig. 5).

Discussion
Our results confirm that the 3S Basin is an important tributary system in the LMB, with impressive biodiversity. For example, the 247 species observed in the 3S Basin comprise 60% of the 411 freshwater fish species found in Cambodia 26 and 22% of the 1148 fish species estimated for the entire Mekong Basin. Moreover, fish communities in the 3S Basin showed marked adaptation to each river basin, characterized by distinct spatial composition and indicator species (Fig. 2, Table 1). Two of the four non-native fish are indicator species, which suggest that the 3S now provides habitat for the non-native species and therefore could threaten native species by competing for resources 24 . Our findings also demonstrate that dams and fish biodiversity are inversely correlated. As the number of dams has increased in the Srepok and Sesan Basins, fish richness, abundance, and diversity have generally declined. On the other hand, fish richness, abundance, and diversity have increased through our study period in the relatively healthy and free-flowing rivers of the Sekong Basin (Fig. 4, Fig. S1).
Average total fish abundance over the study period was not significant different among the three rivers, but migratory fish abundance was significantly higher in the Srepok Basin, while threatened species abundance was significant higher in the Sekong Basin. However, monthly abundance significantly increased through time in the Sekong Basin for total and migratory fish (Fig. S1). Clear evidence of increased fish biodiversity in the Sekong Basin is illustrated by diversity, which significantly increased through time for all fish categories and which accounts for both species richness and their relative abundance 27 (Fig. S1). The opposite results were found for the Srepok and Sesan Basins, specifically the significant decreasing diversity through time for threatened species, and for the total and migratory species, respectively (Fig. S1). www.nature.com/scientificreports/ Our findings highlight the importance of the Sekong Basin as a free-flowing river corridor for fish migration and suggest that the Sekong Basin potentially compensates for environmental changes and habitat loss in the Sesan and Srepok Basins. The Sekong Basin had the highest number of total, migratory, and threatened fish species observed in the 3S Basin (Fig. 3a) and currently provides about 10,330 km of accessible habitat to migratory fish 21 . We counted more total and migratory fish species than past observations 22 . Sekong Basin fish communities included 28 indicator species across the study period, which means that these species were abundant and widespread in the Sekong Basin throughout the study period. The non-native species, Nile tilapia (Oreochromis niloticus), was an indicator species in the Sekong Basin because they were caught only in this basin and were abundant (Table S1). Although the detrimental impacts of this species in the Mekong Basin is not well understood, its presence could lead to the decline of native species 28 . This highlights the need to manage the 3S Basin, especially the Sekong Basin, to prevent future threats to native fish from non-native species. Furthermore, we found six migratory indicator species and two threatened indicator species (the Mekong endemic striped catfish P. hypophthalmus and the Vulnerable M. bocourti) in the Sekong Basin, highlighting the basin's importance to migratory and threatened fish 29,30 .
As dams were built in the Sesan and Srepok Basins, fish diversity decreased. The Sesan is the most dammed river in the 3S Basin (Fig. 1b) and it supported the lowest fish biodiversity, confirming results from a previous study 9 . The Sesan Basin also had the lowest number of migratory, threatened, and indicator species, although in the past it was a migration route for at least 54 migratory species 22 . Five indicator species were recorded in the dry season, one of which is the Mekong endemic Henicorhynchus lobatus, indicating their long-term occurrence with high abundance only in the dry season. The Sesan Basin had no wet season indicator species. However, throughout the study period and regardless of season, the Sesan Basin had only one indicator taxon, Acantopsis Table 1. Indicator species list by river and season in the 3S Basin. Asterisks *P < 0.05, **P < 0.01, ***P < 0.001, m migratory, VU Vulnerable, EN Endangered, CR Critically Endangered, nn non-native, end endemic in Mekong. Values range from 0 to 1, and values near 0 indicate rare or sporadic abundance, while consistently abundant species at many sites near 1. www.nature.com/scientificreports/  www.nature.com/scientificreports/ sp., suggesting that fish species that were present in the beginning of the study period were replaced by different species by the end of the study period. The Srepok Basin was the second-most dammed river in the 3S Basin and supported more fish biodiversity and more indicator species than the Sesan Basin, although richness and diversity for threatened fish species have declined through time. The Srepok Basin historically supported the highest number of migratory species in the 3S Basin 9,21 . However, compared to past records by Baran et al. 22 , the number of migratory fish has declined by 15 species and total species count was reduced by 63 species by 2014. Lower Sesan 2 Dam was completed in 2018 and disconnected the Srepok and Sesan Basins from the Mekong River (LSS2 in Fig. 1). That dam blocked 18,701 km of upstream habitat for fish migration. The presence of Lower Sesan 2 Dam alone underscores the importance of maintaining free-flowing rivers of the Sekong Basin. Our results are the first to demonstrate that fish abundance, richness, and diversity are increasing in the Sekong Basin as other tributaries are disconnected by dams. We have shown that dam construction in the 3S Basin during past decades have altered river ecosystems.
Our results mirror studies in other regions which have shown that dams are an overriding driver of native fish species loss 31 . For example, hydropower development has fragmented river systems in major Andean-Amazonian basins, especially for endemic and migratory species 32 . River fragmentation from future hydropower dam development is anticipated to severely impact fish movement in the Mekong, as well as Amazon, Niger, Congo, and Salween Basins 33 . The combination of climate change, drought intensification, and hydropower dam development is anticipated to harm fisheries and fish biodiversity 34 .
In the Mekong, hydropower is often promoted as a renewable, sustainable source of energy for economic development 14 . However, our study demonstrated that hydropower development causes fisheries and biodiversity loss. This finding contributes to existing research showing that dams degrade ecosystems by altering water flow and quality, blocking fish migration routes, and disturbing fish spawning and feeding grounds 9,35-38 . Free-flowing rivers and intact migration corridors are thus vital to sustaining future food webs, biodiversity, fisheries, and food security in the LMB 39,40 . Protecting migratory fish species is particularly important as they are the food source for threatened birds and freshwater mammals such as the Critically Endangered Irrawaddy Dolphin 41 . Migratory fish also make up about 60% of the catch in the LMB 42 and are central to food security and nutrition because they provide 49-83% of animal protein for people living in the LMB 40 . Moreover, previous research has shown that www.nature.com/scientificreports/ dam construction indirectly leads to deforestation, displaced communities, increased socio-economic conflict, and deterioration of natural resource livelihoods like semi-subsistence farming, wild fishing and forestry [43][44][45][46] . Given the negative impacts of hydropower dams on streamflow and water quality 21,47 , fish biodiversity and abundance, and ecosystem function, it is imperative to include realistic environmental impacts in future hydropower dam planning 21,48 . Our results quantified significant impacts of dams on fisheries in the Sesan and Srepok Basins, and possible compensatory benefits of a largely undammed Sekong Basin. The findings suggest that additional dams in the Sekong Basin, especially the proposed Lower Sekong (LSK) and the Lower Sekong A (LSKA) dams (Fig. 1a), which could be built in the lower mainstem river of Sekong Basin (LSK in Fig. 1a), would likely incur substantial environmental and social costs. Since the Srepok and Sesan Basins have already been dammed, maintaining connectivity in the Sekong Basin to support fish diversity and fisheries is needed and valuable. A largely free-flowing Sekong benefits local communities with diverse freshwater resources and wild fish capture 7,49,50 . Abundant and diverse fish can increase the catch and thus the enough daily protein intake for local people, and provide fisheries-based incomes. Moreover, the free-flowing Sekong can maintain relative natural-seasonal flows that upgrade river integrity and ecosystem services including regulation (e.g. water and erosion regulation, self-purification), provisioning (e.g. water supply and fisheries), and supporting services (e.g. soil formation, nutrient and water cycling). Seasonal agricultural activities can also benefit from the free-flowing www.nature.com/scientificreports/ Sekong because connected rivers ensure seasonal flow availability and support sediment transfer to healthy floodplains 51 , which all are necessary for agricultural production. An alternative to future dams in the Sekong Basin is to promote other alternative renewable resources like solar energy or to re-operate existing dams to maximize hydropower generation without new dams. For example, solar energy combined with centralized power management and trading could enable Thailand, Laos, and Cambodia to meet future energy demand without new hydropower dams 42 . Similarly, operating multiple dams as a system, rather than operating dams independently, could optimize hydropower generation and reduce the need for new dams 52 . In the 3S Basin, 51 dams are already operational, and therefore hydropower re-operation, especially for large dams such as Lower Sesan 2, should be prioritized 53 . Alternative energy sources and optimizing hydropower generation of existing dams will be increasingly important as hydrologic intensification expected with climate change reduces hydropower generation from existing dams and reduces the benefits of proposed dams 54 .

Methods
Study area. The 3S Basin is one of the Mekong's main tributaries and drains southern Laos, Vietnam's central highland, and northeastern Cambodia, joining the Mekong mainstem in Stung Treng Province, Cambodia (Fig. 1a). The 3S Basin is important to Mekong River hydrology, in particular Tonle Sap's flood pulse and seasonal reverse flow, because it contributes up to 25% of the Mekong's annual discharge and sediment load 55 . Like the LMB, the hydrology of the 3S Basin has well-defined wet (May-October) and dry (November-April) seasons. More than 80% of annual rainfall occurs in the wet season. Across the 3S Basin, mean annual rainfall varies greatly from ~ 1500 mm in the lower reaches to ~ 2500 mm in the southern headwaters and > 3000 mm in the northern headwaters 56 .
Data collection. Fish monitoring was conducted from June 2007 to May 2014 in the 3S Basin. Each of the three monitoring sites, that could extend a few kilometers in length, was located in each 3S river adjacent to the village where participating fishers were based (Fig. 1a). The three sampling sites remained relatively unchanged over the course of study. At each sampling site, three professional fishers, who are 'typical' medium-scale fulltime fishers, were trained to record their daily catches by responsible researchers from the Inland Fisheries Research and Development Institute (IFReDI), Cambodia Fisheries Administration. After training, the fishers were tested to ensure their reliable ability in setting up fishing gears, identifying, measuring, recording and taking photos of fish species. For more detailed fisher training and testing, we refer to Ngor et al. 57 .
To reduce the effect of sampling effort among study sites, only fish catches from stationary gillnets (120 ± 50 m long, 2-3.5 m high, with mesh sizes of 3-12 cm) each day (12 ± 2 h per day) were used in this study. Fish were identified and counted using a fish species list (∼ 900 species) from the MRC Mekong Fish Database, FishBase, and other existing literature 1,58 . A fish photo book containing 272 common Mekong's fishes including 20 nonnative species with species code, local and scientific names was also provided to each participating fishers to aid fish identification and fish recording on the data sheet in the field 25 . The fish sampling procedures used in this study follows the MRC standard sampling procedures for Fish Abundance and Diversity Monitoring in the LMB described in Ngor et al. 57 .
Data analysis. In this study, daily fish samples were averaged across all fishers for monitoring sites to daily mean samples, where were then aggregated into monthly fish community data. Monthly data were used to capture fish biodiversity response to seasonal hydrology and climate. Daily data added noise, but did not change monthly trends. Descriptive statistics including mean, standard deviation, minimum and maximum monthly fish species richness and abundance (total individuals catch of all species), and differences between the monitoring sites are provided in Table 2.
To assess spatial-seasonal changes in fish biodiversity in the 3S Basin, a principal component analysis (PCA) was performed using the rda() function of the "vegan" R package 59 . We used Hellinger-transformed data of the total monthly fish community using the decostand() function of the "vegan" R package 59 . This transformation reduced the impact of the highest fish abundance values 60 . PCA reduces dimensionality and summarizes the main features of a dataset. For example, data points may separate into clusters when samples are taken from different seasons or ecosystems or gradients of sites when a given biotic community is adapted to different conditions of the ecosystems 60 .
We also determined indicator fish species for each 3S river. An indicator fish species here is a species that occurs consistently through time, and more abundant compared to other species. Indicator species value ranges from 0 to 1. Indicator values near 0 indicate rare or sporadic abundance. When a species is consistently abundant in sites and is found at all sites of that group, and therefore has a large mean abundance within the group, its indicator value is 1. Smaller indicator values suggest that a species is not an indicator species; however, smaller indicator values do not necessarily preclude species from being indicators because the value also depends on the abundance of the species caught in each sample. We used the "indicspecies" R package 61 to determine the indicator fish species and compute and test their indicator values against a P value of 0.05 based on a permutation test 60 . If any species has a significant indicator value, it is considered an indicator species. An ecosystem that supports a high number of indicator species can provide information on the characteristics and environmental conditions of the habitats and sites that they share 61 . A species was considered an indicator when its occurrence and mean abundance in the study site remained high throughout the study period, implying favorable ecological conditions at that site throughout the study period 60 . In this regard, a river that supports a number of fish indicator species may indicate that the river's ecological conditions are favorable and stable over time. Otherwise, the fish species www.nature.com/scientificreports/ occurring in the river are merely visitors or do not occur consistently in the river due to hydrological or water quality fluctuations, habitat alteration, or river fragmentation. Fish biodiversity was grouped into three fish categories: (1) total species, (2) migratory species and (3) threatened species. Each category differs in sensitivity to river conditions, river connectivity, and imminent risk to biodiversity. Migratory fishes undergo regular, predictable movements from one region to another within a larger system (www. fishb ase. se), and threatened fish species are listed as Critically Endangered (CR), Endangered (EN) and Vulnerable (VU) 62 . The fish biodiversity metrics such as species richness, abundance and diversity (refereeing to Shannon's diversity index) were used for all three fish categories for each 3S river. Species richness and species abundance refer to the number of different species and the total individuals of all species per sample or per given area, respectively. The Shannon's diversity index is computed following the equation: diversity = -1Σ [p i × ln (p i )] 27 . Here, p i is the proportion of individuals of the ith species in a given community, and is computed as p i = n/N, where n is the individuals of the ith species and N is the total number of individuals in the whole community. The metrics were compared to understand whether fish communities in the Sekong, Srepok, and Sesan Basins are distinct from one another. The comparison was based on one-way ANOVA using the aov() function in the "stats" R package, or Kruskal-Wallis tests using the kruskal.test() functions 63 . The former was applied when residuals of the Shapiro-Wilk test were normal; otherwise, the non-parametric test was used. When significant differences were detected, we performed Wilcoxon tests using the wilcox.test() function of the "stats" R package 63 to evaluate the differences in fish biodiversity metrics among river basins.
To investigate the temporal trends in fish biodiversity, we employed linear models to regress the monthly biodiversity metrics of the three fish categories against time using the lm() function of the "stats" R package 63 . In total, there were 84 monthly data points of fish biodiversity. To quantify the impacts of dam development, each metric of the three fish categories from each river was averaged to yearly data and then regressed against the yearly number of hydropower dams that were constructed in the 3S Basin. This produced eight data points for the 2007-2014 sampling period. We used yearly data because data for the month that hydropower dams became operational in unavailable. All analyses were conducted using the R statistical programing language 63 and values of P < 0.05 were considered significant. The map in this study was created using ArcMap version 10.4.1. All figures were generated using functions in the "stats" R package 63 and the "ggplot2" R package 64 .
Ethical statements. All field and measurement protocols were approved by the Inland Fisheries Research and Development Institute (IFReDI), Cambodia Fisheries Administration, and all methods were carried out in accordance with relevant guidelines and regulations of IFReDI. Fish catch and handling after all scientific measurement and records was granted to local fisher for consumption and sale by Fisheries Capture Committee, and the need for ethics approval is deemed unnecessary according to the regulations of IFReDI.